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The Fermi Gamma-ray Burst Monitor (GBM) detects gamma-rays in the energy range 8 keV - 40 MeV. We 
developed a new background fitting process of these data, based on the motion of the satellite. Here we 
summarize this method, called Direction Dependent Background Fitting (DDBF), regarding the GBM triggered 
catalog. We also give some preliminary results and compare the duration parameters with the 2-years Fermi 
Burst Catalog. 



1. Introduction 



one can plot a GBM lightcurve as shown in Fig. [TJ 



Fermi has specific proper motion when surveying 
the sky. It is designed to catch gamma-ray bursts in 
an effective way. However, bursts can have a varying 
background, especially in the case of the Autonomous 
Repoint Request (ARR). Modeling this with a poly- 
nomial function of time is not efficient in many cases. 
Here we present the effect of these special moving fea- 
ture, and we define direction dependent underlying 
variables. We use them to fit a general multidimen- 
sional linear function for the background. 

Note that here we give a short summary of the 
method and some preliminary results. The detailed 
description of the Direction Dependent Background 
Fitting (DDBF) algorithm and the final results will 
be pubhshed soon [Szecsi et al.||2013 . 



3. Direction Dependent Bacltground 
Fitting (DDBF) 

Analyzing ancillary spacecraft and other directional 
data we have found the following (x^) variables, which 
seem to contribute to the variation of the background: 
celestial distance between burst and detector orienta- 
tion, celestial distance between Sun and detector ori- 
entation, rate of the Earth-uncovered sky and time 
Szecsi et aL]|2013 



We use the method of General Least Square for 
multidimensional fitting of the yi counts to the cor- 
responding (x) explanatory variables. The maximum 
likelihood estimate of the model parameters is 
obtained by minimizing the quantity of 



2. Fermi lightcurves with varying 
bacltground 

Fermi's slewing algorithm is quite complex. It is 
designed to optimize the observation of the Gamma- 
Ray Sky. In Sky Survey Mode, the satellite rocks 
around the zenith within ±50° , and the pointing alter- 
nates between the northern and southern hemispheres 
each orbit |Fitzpatrick]|2011[ [Meegan et alT]|2009| . In 
Autonomous Repoint Request (ARR) mode, it turns 
toward the burst and stays there for hours. 12 Nal 
detectors are placed such a way that the entire hemi- 
sphere is observable with them at the same time. 

The GBM data, which we use in our analysis (called 
CTIME), are available at 8 energy channels, with 
0.064-second for triggered and 0.264-second resolution 
for non-triggered mode. The position data is avail- 
able in 30-second resolution. This data were evenly 
proportioned to 0.256-second and 0.064-second bins 
using linear interpolation, in order to correspond to 
the CTIME data [Meegan et al. 2009] . 

Using GBM data, 1-second bins and summarizing 
the counts in the channels between 11.50-982.23 keV, 



X — ^i=i 



Hi - S^f^^afeXfc(xj) 



The matrix of A and vector b are: 



CT,; cr,; 



(1) 



(2) 



Minimizing leads us to the following equation: 



a= (A'^A)-iA'^b 



(3) 



, where A"^ means the transpose of A, and the ex- 
pression (A-^A)~^A-^ are called generalized inverse 
or pseudoinverse of A. For calculating the pseudoin- 
verse of the design matrix A we used Singular Value 
Decomposition (SVD): A = USF^. Using U and V, 
the pseudoinverse of A can be obtained as 



pinv{A) = {A'Ay^A^ = VS^'U 



3-1ttT 



(4) 



Since we want to have a method for all the Fermi 
bursts, we define our model to be quite comprehen- 
sive. Let us have y(xi) as the function of x^ = 
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Figure 1: Lightcurve of the Fermi burst 091030613 measured by the 3rd GBM-detector, without any background 
filtering, with 1-second bins. The grey line is a fitted polynomial function of time of order 3. 



{xf'\x^j^\xf'\x^f^) of order 3, so the basis functions 
Xfc(xi) (and columns of the design matrix) consist of 

every possible products of the components ^ up to 
order 3. That means that we have M = kmax = 35 
basis functions and oi, 02-. .035 free parameters. It is 
sure that we do not need so many free parameters to 
describe a simple background. Computing the pseu- 
doinverse we need the reciprocal of the singular values 
in the diagonals of S~^, but if we compute the pseu- 
doinverse of A, the reciprocals of the tiny and not im- 
portant singular values will be unreasonably huge and 
enhance the numerical roundoff errors as well. This 
problem can be solved defining a limit value, below 
which reciprocals of singular values are set to zero 
Szecsi et ar]|2013| 



One cornerstone of the fitting algorithm described 
above is the definition of the boundaries which divide 
the interval of the burst and the intervals of the back- 
ground. In this work, we follow the common method 
of using user-selected time intervals Paciesas et al. 



2012 . But unlike Paciesas et al. 2012 , using the 



position data allows us to fit the total CTIME file in- 
stead of selecting two small intervals around the burst 
[Szecsi et al.||2013| . 



4. Model selection 



ues from the model are normally and independently 
distributed, every model has a value AlCfc so that 



AlCfe = TV • log : 



RSSh 



N 



2-k 



(5) 



, where RSSk is the Residual Sum of Squares from the 
estimated model {RSS = (y; — y{xi,k))'^), N is 
the sample size and k is the number of free parameters 
to be estimated. Given any two estimated models, the 
model with the lower value of AIC^ is the one to be 
preferred. Given many, the one with lowest AIC^ will 
be the best choice: it has as many free parameters as 
it has to have, but not more. 

We loop over the pseudoinverse operation and 
choose Skk as the limit of singular values in the fcth 
step, and compute the corresponding AIC^. The num- 
ber of singular values which minimalize the AlCfc as a 
function of k will be the best choice when calculating 
the pseudoinverse .Szecsi "eraL][20T3l . 

As an example we analyze the lightcurve of 
GRB 090113778. GRB 090113 is a long burst with 
T§g*=17.408±3.238 s in the GBM Catalogue, here we 
show Detector data: 

Fig. [2] has some extra counts around 400 and 600 
seconds. Both of them can be explained with the vari- 
ation of the underlying variables, that is, the motion 
of the satellite. These are not GRB signals! 



The Akaike Information Criterion (AIC) is a com- 
monly used method of choosing the right model to 
the data [Akaike 1974 . Assume that we have M 
models so that the fcth model has fc free parameters 
(fc = 1...M). When the deviations of the observed val- 



5. Error analysis 

The DDBF method is too complicated to give a 
simple expression for the error of Tqq using general 
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Figure 2: Up: Lightcurve of the Fermi GRB 090113778 
measured by the triggered GBM-detector '8', and the 
fitted background with a grey line. Burst interval (sees): 
[-20:40]. Down: Underlying variables (absolute values). 



Figure 3: Down left: Akaike Information Criterion. A 
model with 12 singular values was selected. Down light: 
Cumulated lightcurve of GRB 090113778. 



5.2. GRB 091030613 



rules of error propagation. We therefore repeated the 
process for 1000 MC simulated data. 



5.1. GRB 090113778 



Distribution of the Poisson-modified Tgg and T50 
values are shown in Fig. g] for GRB090113778. 

Fig. [4] shows the Monte Carlo simulated distribu- 
ton of the duration values of GRB 090113778. Based 
on this, we give confidence intervals corresponding to 
68% (Icr level). The result is given in Table [l] Table |l) 
also contains some of our other, preliminary results 
using DDBF. For more details, see our forthcoming 
paper Szecsi et al. 2013 . 



Distribution of the Poisson-modified Tgo and T50 
values are shown in Fig. [5] for GRB091030613. 

Fig. [4] shows two significant peaks around 22 and 
47 seconds. The first peak at 22 seconds corresponds 
to the measured Tgo value. However, in some cases of 
the Poison noise simulation, the measured Tgo value is 
systematically longer: that is because this burst has 
a little pulse around 47 seconds (see Fig. [I]). There is 
no sign of this second peak in the Tso distribution, as 
it is more robust. 



5.3. Preliminary results 



In Figures [6p0] show the MC simulated Tgo and T50 
values for the GRBs given in Table |l] 
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Burst 


79()(s) 


errors (s) 


rpcatalog ! \ 
-'90 \°l 


T5o(s) 


errors (s) 


rpcatalog / \ 

-'so y^) 


081009690 


176.228 


+ 1.357 


—9.477 


176.191 


15.852 


+3.006 


—2.350 


25.088 


090102122 


29.756 


+2.971 


-1.198 


26.624 


10.859 


+0.531 


-0.556 


9.728 


090113778 


19.679 


+ 10.883 


-6.421 


17.408 


6.408 


+0.498 


-0.344 


6.141 


090618353 


103.338 


+3.842 


-6.725 


112.386 


22.827 


+2.201 


-1.530 


23.808 


090828099 


63.608 


+1.467 


-1.652 


68.417 


11.100 


+0.198 


-0.194 


10.752 


091030613 


22.609 


+13.518 


-4.522 


19.200 


10.770 


+0.388 


-0.424 


9.472 


100130777 


80.031 


+3.755 


-3.485 


86.018 


32.340 


+0.931 


-1.363 


34.049 



Table I Preliminary results and confidence intervals. aj^j rpcatalog from the GBM Catalogue 
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Figure 5: Distibution of the Tgo (left) and Tso (right) 
values obtained from the MC simulated data for Fermi 
burst 091030613. 
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Figure 6: Distibution of the Tgo (left) and T50 (right) 
values obtained from the MC simulated data for Fermi 
burst 081009360, GBM detector '8'. 



Figure 4: Distibution of the T90 (up) and Tso (down) 
values obtained from the MC simulated data for Fermi 
burst 090113778. 



6. Summary 

We summarized the Direction Dependent Back- 
ground Fitting (DDBF) algorithm which is designed 
to filter the background of the Fermi lightcurves 
on a longer timescale (2000 seconds of the CTIME 
datafile). This technique is based on the motion 
and orientation of the satellite. The DDBF con- 
siders the position of the burst, the Sun and the 








Figure 7: Distibution of the Tgo (left) and T50 (right) 
values obtained from the MC simulated data for Fermi 
burst GRB 090102122, GBM detector 'a'. 
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Figure 8: Distibution of the Tgo (left) and Tso (right) 
values obtained from the MC simulated data for Fermi 
burst 090618353, GBM detector 7'. 



we computed three physically meaningful underlying 
variables, and fitted a 4 dimensional hypersurface on 
the background. Singular value decomposition and 
Akaike information criterion was used to reduce the 
number of free parameters. 

One of the main advantage of the DDBF method 
is that it considers only variables with physical mean- 
ing. Furthermore, this method can fit the total 2000 
sec CTIME data as opposed to the currently used 
methods. This features are necessary when analyzing 
long GRBs and precursors, where motion effects influ- 
ence the background rate sometimes in a very extreme 
way. Therefore, not only Sky Survey, but ARR mode 
GRB's can be analyzed, and a possible long emission 
can be detected. 

A detailed description of the DDBF method, the 
first resulst, the comparison with the currently used 
methods, and analysis of the ARR cases will be pub- 
lished soon Szecsi et"ar]|2013 . 



Figure 9: Distibution of the Tgo (left) and Tso (right) 
values obtained from the MC simulated data for Fermi 
burst 090828099, GBM detector '5'. 



Acknowledgments 

This work was supported by OTKA grant K077795, 
by OTKA/NKTH A08-77719 and A08-77815 grants 
(Z.B.). The authors are grateful to Aron Szabo, Peter 
Veres for the valuable discussions. 




Figure 10: Distibution of the Tgo (left) and Tso (right) 
values obtained from the MC simulated data for Fermi 
burst 100130777, GBM detector 'b'. 



Earth as well. Based on these position information, 
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